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Abstract 

We constrain the true spin-orbit alignment of the KOI-89 system by numerically fitting the two Kepler 
photometric lightcurves produced by transiting planets KOI-89.01 and KOI-89.02. The two planets have 
periods of 84.69 days and 207.58 days, respectively. We find that the two bodies are low-density giant 
planets with radii 0.45 ± 0.03 Rjup and 0.43 ± 0.05 Rjup and spin-orbit misalignments 72° ± 3° and 73°^5^, 
respectively. Via dynamic stability tests we demonstrate the general trend of higher system stability with 
the two planets close to mutual alignment and estimate their coalignment angle to 20° ± 20° - i.e. the 
planets are misaligned with the star but may be aligned with each other. From these results, we limit KOI- 
89’s misalignment mechanisms to star-disk-binary interactions, disk warping via planet-disk interactions, 
planet-planet scattering, Kozai resonance, or internal gravity waves. 

Keywords: techniques: photometric — planets and satellites: fundamental parameters — planets and satel¬ 
lites: formation — planets and satellites: physical evolution — planets and satellites: individual: 
KOI-89.01, KOI-89.02 — stars: rotation 


1. INTRODUCTION 


Recent studies show that exoplanetary systems around 
early-type stars display a high diversity in their fundamen¬ 
tal characteristics with varied planet radii, planet densities, 
periods, eccentricities, and inclinations. (Brandt et al. 2014; 
Borucki 2012; Howard 2013). Exoplanet orbits have notably 
caried spin-orbit alignments with highly misaligned and even 
retrograde orbit geometries (Bourrier & Hebrard 2014; Al¬ 
brecht et al. 2012; Winn et al. 2011). At the time of this 
work, most known misaligned systems are short-period; only 
HD80806b (Naef et al. 2001), Upsilon Andromedae (Deitrick 
et al. 2015), and Kepler 56 (Huber et al. 2013) have peri- 
Dds > 10 days.^ This work adds to the list of long-period 
spin-orbit misaligned planets KOI-89.01 and KOI-89.02. 

There are several methods for determining various aspects of 
a system’s alignment, including gravity-darkening (J. Barnes 
et al. 2011; Ahlers et al. 2014), the Rossiter-McLaughlin ef¬ 
fect (Ohta et al. 2005), Doppler tomography (Gandolfi et al. 
2012), asteroseismology (Ghaplin et al. 2013; Van Eylen 
et al. 2014), photometric amplitude distribution (Mazeh 
2015), and stroboscopic starspots (Desert et al. 2011; Hi- 
rano et al. 2012). We applied the gravity-darkening method 
first suggested by J. Barnes (2009) and later applied to Ke¬ 
pler Object of Interest (KOI) 13 (J. Barnes et al. 2011) and 
KOT2138 (J. Barnes et al. 2015). This method constrains 
both the star’s polar tilt toward/away from the observer 
(stellar obliquity), and the planet’s misalignment angle as 
seen relative to the observer (sky-projected alignment). 

The gravity-darkening effect, first predicted by Von Zeipel 
(1924), results in a pole-to-equator gradient in stellar lumi¬ 
nosity driven by rotation. As an object transits a gravity- 
darkened star, it can move across areas of unequal bright¬ 
ness; this luminosity gradient can affect the lightcurve in 

^ http://www.physics.mcmaster.ca/'rheller/ provides a list of spin- 
orbit misaligned planets. 


various ways, depending on its transit geometry (J. Barnes 
2009). J. Barnes et al. (2011), J. Barnes et al. (2013), Zhou 
& Huang (2013), and Ahlers et al. (2014) all showed that 
the asymmetry in such lightcurves (or lack thereof) can be 
utilized to constrain the spin-orbit alignment of a transiting 
system. 

The causes of frequent misalignment around fast rotators are 
still under investigation. The underlying issue is that planets 
probably do not form with initially misaligned orbits - their 
angular momenta must be conserved with the stellar nursery 
they formed in. In this case, such planets must migrate to 
their misaligned positions. KOI-89 is one such system that 
does not conform to the traditional nebular hypothesis. 

There are several ideas for processes that might create spin- 
orbit misalignment. Lai et al. (2011) and Spalding & Baty¬ 
gin (2014) demonstrated that magnetic torques can push the 
stellar spin axis away from the circumstellar disk’s angular 
momentum vector over very long timescales. This would 
specifically explain spin-orbit migration in very young sys¬ 
tems with late-type stars, where stellar magnetic fields are 
strongest. Rogers et al. (2012) showed that internal gravity 
waves can produce angular momentum transport between 
the convective interior and radiative exterior of early-type 
stars that turn the stellar spin axis away from the system’s 
invariant plane. 

There are also several ideas that explain how spin-orbit mi¬ 
gration might develop via more dynamic means. Libert & 
Tsiganis (2009) discussed Kozai resonance in a 2-planet sys¬ 
tem and its effects on mutual inclination. This is almost 
certainly the origin of misalignment for HD 80806b (Naef 
et al. 2001). Ghatterjee et al. (2008), Ford et al. (2005), Ray¬ 
mond et al. (2008), and Nagasawa et al. (2008) all demon¬ 
strated how planet-planet scattering can drive misalignment 
in a multiplanet system. Levison et al. (1998) showed that 
planet-embryo collisions during planet formation can lead 
to high mutual inclination. Terquem (2013), Teyssandier 
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et al. (2013), and Batygin (2012) analyze gravitational disk¬ 
warping events that lead to misalignment. 

Winn et al. (2010) showed a correlation between hot stars 
(2eff > 6250) and misalignment. Batygin & Adams (2013) 
showed an interdependence between stellar rotation rates 
and spin-orbit misalignments. These works imply that a 
large number of planets orbiting early-type stars are com¬ 
monly misaligned. 

Huber et al. (2013) employed asteroseismology to measure 
the the stellar obliquity of multiplanet system Kepler 56, 
and showed that spin-orbit misalignment is possible in mul¬ 
tiplanet systems with low-mass, long-period planets. Beno- 
mar et al. (2014) found mild misalignment in Kepler-25 via 
a joint analysis of asteroseismology, lightcurve analysis, and 
the Rossiter-Mclaughlin effect. 

This work provides another example of a long-period multi¬ 
planet system with significant misalignment: KOT89. In §2, 
we outline our data preparation process and list previously 
known system characteristics. In §3, we introduce new tech¬ 
niques to the J. Barnes et al. (2011) fitting method. In §4, 
we show our best-fit parameters and constraints on misalign¬ 
ment. We test KOI-89’s dynamic stability and constrain the 
coalignment angle between the two orbits in §5. In §6, we 
discuss possible formation and migration mechanisms for the 
KOT89 system, as well as test the dynamic stability of the 
system in order to constrain the planets’ mutual alignment. 

2. OBSERVATIONS 
2.1. Data Preparation 

The Mikulski Archive for Space Telescopes (MAST) Kepler 
Input Catalog (KIC) provided the Kepler photometry that 
we analyze for the KOI-89 system. We employ each of the 16 
available quarters of KIC data, combining them into a single 
dataset. We choose to only incorporate long cadence data 
(30-minute integrations) because both planets have transit 
durations of over 12 hours. Therefore ingress and egress are 
well-sampled by the 30 minute time cadence, and inclusion of 
short (1-minute) cadence data would not provide additional 
constraints. 

After concatenating all available long cadence photometry, 
we apply a median box filter of 44 hours (three times KOT 
89.01’s transit duration) to reduce long-term astrophysical 
and instrumental variability. Figure 1 displays the filtered 
time-series. We then identify which transits correspond to 
which transiting body based on their KIC orbital periods, 
and separate them accordingly into individual datasets. 

We adjust the center-of-transit times of each transit 
lightcurve according to their measured transit timing vari¬ 
ations (TTV) (Rowe et al. 2014). We perform this adjust¬ 
ment for each measured TTV in Rowe et al. (2014), includ¬ 
ing thirteen transits for KOI-89.01 and five transits for KOT 
89.02. These transits exclude the double transit identified 
in Figure 1. 

With all TTV accounted for and the individual transits 
evenly separated by 84.69 days and 207.58 days, respec¬ 
tively, we fold all KOT89.01 transits on top of the epoch 
34960800 ± 400 s transit and fold all KOT89.02 transits on 
top of the epoch 25041400 ± 700 s transit. We then combine 
the two resulting lightcurves back in a single dataset. With 
the Kepler photometry represented by a single lightcurve 
with two transit events, we bin the data at 15 minutes to 
improve the computation time of our fit. We determine the 
error bars of the binned data from the standard deviation 
of the flux values in the bin. 
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Figure 1. Kepler photometry of the KOT89 system. The verti¬ 
cal length of the data points represent their uncertainties. Transits 
of two planets are visible, displaying periodicities of 84.69 days and 
207.58 days, respectively. In mid-2011 (arrow), both planets transited 
simultaneously. 

2.2. Previously Measured KOT89 Parameters 

We build our work upon previous research of the KOT89 
system. We obtain the KOT89 stellar mass, stellar temper¬ 
ature, and transit periods from the Community Follow-up 
Observation Program (CFOP). Figure 2 shows the spectro¬ 
scopic determination of '!;sin(t). We list these and other 
relevant system parameters in Table 1. 


Parameter 

Values 

Pi 

84.69 days 

P2 

207.58 days 

PeB 

7717 ± 225 K 

M* 

1.965 ±0.256 M0 

V sin(2) 

90 km/s 

Vmag 

11.731 

KOI-89.01 SNR 

93.8 

KOI-89.02 SNR 

68.1 


Table 1 

Previously measured parameters of the KOI-89 system. We 
incorporated all parameters as assumed values when fitting the 
KOI-89 lightcurve. 

Rowe et al. (2014) confirmed 715 new systems - including 
KOT89 - via multiplicity. The two planets have a period 
ratio near the 5:2 mean-motion resonance (2.45). Follow¬ 
up observations of these phenomena could confirm/deny the 
existence of additional orbiting bodies, and could further 
constrain this system’s formation and evolution. 

3. MODEL 

We update the transitfitter program (J. Barnes 2009) 
to fit multiple-planet transiting systems. The Levenburg- 
Marqhardt minimization technique remains the fitting 
method, but now treinsitf itter can constrain the orbital 
parameters of additional transiting bodies. These extra pa¬ 
rameters cause an increased sensitivity to the photomet¬ 
ric signal-to-noise ratio, so minimization must be ap¬ 
proached with additional caution. 

The individual parameters of additional planets have the 
same degeneracies as a single-planet fit. There is a degen- 
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K00089, obs = TR, 2013-11 -15_03h07tn41 s, hjd = 2456611.634598, ap = 23, expt = 720 s, SNRe = 40.4, ver = 2.5, SPCv = 2.2 



Figure 2. Spectroscopic determination of Dsin( 2 ) for the KOI-89 system, measured with the Tillinghast Reflector Echelle Spectrograph on the 
1.5 m telescope at the Whipple observatory. These data were provided by David Latham of the Kepler Science Team and retrieved from the 
Community Follow-up Observing Program (CFOP). 

Ip, producing a coherent physical system. In our fitting 
model, Ip is particularly sensitive to the lightcurve geome¬ 
try; using multiple lightcurves simultaneously allows ip to be 
constrained by more data points, decreasing its uncertainty. 

4. RESULTS 
4.1. Lightcurve Fits 

The lightcurves of KOI-89.01 and KOI-89.02 (Figure 4) dis¬ 
play unusual shapes. KOI-89.01 has the asymmetry ex¬ 
pected of a misaligned body orbiting a fast-rotator (J. 
Barnes 2009) . KOI-89.02 does not display this asymme¬ 
try, possibly due to lower photometric precision. Both tran¬ 
sits show sloped ingresses/egresses and entirely non-constant 
transit depths, producing dominant V-shaped lightcurves. 
A typical lightcurve is symmetric with a steep ingress and 
egress with a relatively flat bottom, rounded only by limb 
darkening. 

KOI-89’s V-shaped lightcurves can arise in one of two 
ways. First, planets only grazing their star during tran¬ 
sit rather than fully eclipsing it block constantly chang¬ 
ing sky-projected areas. This effect creates sloped in¬ 
gresses/egresses. However, this situation is improbable for 
KOI-89 as both planets would require similar, high impact 
parameter values despite having significantly different semi¬ 
major axes. 

The second way KOI-89 could generate V-shaped lightcurve 
geometries is by having a gravity-darkened star with a very 
high stellar obliquity ip - i.e. pole-on. In this case, the plan¬ 
ets transit near a stellar pole and the gravity-darkened equa¬ 
tor surrounds the outer edge of the star. The limb-darkening 
and gravity-darkening effects combine together to create a 
significant center-to-edge luminosity gradient. At ingress, 
the planet blocks a continuously increasing total flux as it 
moves closer toward the center of the star, and vice versa 
during egress. This produces a V-shaped transit lightcurve 
for each planet (J. Barnes (2009), Figure 4), consistent with 
the lack of the typical ingress-egress asymmetry expected in 
a misaligned gravity-darkened transit. We test for the pos¬ 
sibility of grazing transits by fitting the system with impact 
parameters nearly at and slightly above 1.0. We find that we 
can not match the system’s lightcurve with grazing transits: 
such an event can not reproduce the proper ingress-egress 
asymmetry seen in KOI-89.01. 


eracy between eccentricity and stellar radius, for instance: 
an eccentric orbit can have the same transit duration as a 
circular orbit around a smaller star. Also, fitting transit 
lightcurves in isolation cannot determine stellar mass, so we 
apply an assumed stellar mass from CFOP and fit the ec¬ 
centricity around it. A transit around a fast-rotator has 
degenerate limb-darkening and gravity-darkening effects in 
the case of high stellar obliquity. We discuss this degeneracy 
further in §4. 



Figure 3. Top left: The stellar obliquity is defined as the axial 
tilt toward/away the plane of the sky. Top right: the sky-projected 
alignment A is the misalignment angle seen from Kepler’s point of 
view. Bottom left: the spin-orbit alignment tp is the angle between the 
plane of an orbit and the star’s equatorial plane. Bottom right: the 
coalignment angle a is the angle between two orbit planes. 

Fitting for both planets simultaneously rather than fitting 
each lightcurve separately forces the stellar radius and 
the stellar obliquity ip (shown in Figure 3) to agree for both 
planets, which has two advantages. First, a simultaneous 
two-planet fit requires two fewer parameters to be fit, in¬ 
creasing the overall accuracy and decreasing computation 
time of the other parameters. Second, it applies lightcurve 
geometries of both planets toward the constraint of i?* and 
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Figure 4 shows our best-fit lightcurve using grazing transits 
in blue. We apply grazing transits to both spherical and 
gravity-darkened models. We hold the stellar obliquity at 
zero in the gravity-darkened model to test the system for 
possible spin-orbit alignment. The poor fit of Xreduced = 
1.94 (adjusted to account for holding the stellar obliquity 
constant) motivates us to investigate a model with a high 
stellar obliquity and rapid stellar rotation. 

Using the Levenberg-Marqhardt minimization technique, 
we fit for thirteen parameters: 

• The stellar equatorial radius (i?*) 

• The stellar obliquity (tp) 

• The stellar normalized flux (Fg) 

• The radii of KOI-89.01 and KOI-89.02 (Rp,, Rp^) 

• The inclinations of KOI-89.01 and KOI-89.02 (U, 12) 

• The sky-projected alignments (Ai, A2) 

• The two orbits’ eccentricities (ei , 62) 

• The center-of-transit times (Tq^, Tq^) 

We display the best-fit lightcurve of our gravity-darkened 
model in Figure 4 as the red line. 


Parameter 

Best Fit Values 

^reduced 

1.52 

R* 

2.3 ±0.2 Rq 

ifj 

69° ± 3° 

Cl (fixed) 

0.56 

C 2 (fixed) 

-0.15 

/3 (fixed) 

0.25 

Po 

1.000009±3* 10-® 

Prot (derived) 

S.Sllj;® hr 

/* (derived) 

n iq+0.04 
'4'ry_o.o3 

Ppi 

0.45 ± 0.03 Rjup 

Pp 2 

0.43 ± 0.05 Rjup 

ei > 

0.056 ±0.019 

62 > 

0.50 ±0.09 

h 

89.340° ± 0.05° 

n 

90.64° ± 0.06° 

bi (derived) 

^■°4_o.07 

62 (derived) 

y-O'-o.ir 

n. 

34960800 ± 400 s 

P 02 

25041400 ± 700 s 

Al 

-32° ± 11° 

A 2 

-32° ± 40° 

ipi (derived) 

72° ± 3° 

ip 2 (derived) 

.j. 3 °+ii 


Table 2 

Best-fit results for the KOI-89 system. We calculated stellar period of 
rotation Prot from usin( 2 ), M*, and -i/;. We derived the stellar 
oblateness / from the Darwin-Radau relation. The impact 
parameters bi and 62 were found using Pi and P 2 , ii and 22 , and R*. 
We set our limb-darkening parameters ci = ui U 2 and C 2 = ui — U 2 
according to Sing (2010). 


4.2. Testing For TTV Systematics 

We test the TTV ephemeris reported in Rowe et al. (2014) 
for systematic errors by fitting the KOI-89 series as two 
epochs. The first epoch is comprised of KOI-89.01’s first 
seven transits and KOI-89.02’s first two transits. The sec¬ 
ond epoch is comprised of KOI-89.Ol’s remaining six transits 
and KOI-89.02’s remaining three transits. For each half, we 


adjust all transits with respect to their TTV and fold the 
transits in the same fashion as described in §2.1. 

We apply our gravity-darkening model to both epochs and 
find that the resulting parameters of each fitted dataset have 
overlapping Icr values with our best-fit values using the full 
timeseries (Table 2). We therefore detect no evidence of sys¬ 
tematics in the TTV ephemeris listed in Rowe et al. (2014). 

4.3. Testing Limb Darkening Assumptions 

Limb-darkening has traditionally been problematic for the 
the gravity-darkening technique. Via Doppler Tomography, 
Johnson et al. (2014) measured the sky-projected alignment 
of KOI-13.01 and found it to differ significantly from the 
gravity-darkening measurement performed in J. Barnes et al. 
(2011). Masuda (2015) proposed a solution to this discrep¬ 
ancy by demonstrating that the gravity-darkening model 
produces concurring measurements with Doppler Tomogra¬ 
phy when using a nonzero second quadratic limb-darkening 
term (02). This limb-darkening term is also a possible expla¬ 
nation for KOI-368’s different spin-orbit misalignment val¬ 
ues measured in Zhou & Huang (2013) Ahlers et al. (2014). 
These works motivated us to update our gravity-darkening 
model to include both quadratic terms, ci and C2. 

KOI-89’s very high stellar obliquity brings about an addi¬ 
tional challenge in resolving limb-darkening. With the stel¬ 
lar pole near the center of the sky-projected stellar disk, the 
gravity-darkening and limb-darkening luminosity gradients 
behave nearly identically and are essentially additive. The 
resulting combined effects on a transit lightcurve are there¬ 
fore degenerate with a stellar obliquity near 90°. 

KOI-89’s spectroscopically-determined effective tempera¬ 
ture of 7717±225 K corresponds to an approximate range of 
0.55 to 0.57 for ci and —0.165 to —0.135 for C2 (Sing 2010). 
We test the robustness of our assumed limb-darkening val¬ 
ues by rehtting using (0.55,-0.135) and (0.57,-0.165) for 
Cl and C2, respectively. 

Applying the assumed (ci, C2) values (0.55, —0.135), we mea¬ 
sure a slight increase in KOI-89.02’s impact parameter; how¬ 
ever, this increase vanishes when adjusting the gravity- 
darkening value P to match Affair’s value of 0.19 (Mon- 
nier et al. 2007). We detect no significant changes in our 
best-fit results when employing the limb-darkening values 
(0.57, —0.165). We cannot resolve the accuracy of our limb- 
darkening parameters without higher precision data, and 
therefore elect to apply the assumed values of /3 = 0.25 
(Von Zeipel 1924) and (ci,C2) values of (0.56,-0.15) (Sing 
2010 ). 

4.4. Eccentricities 

We constrain the lower limits of eccentricity to 0.056 ± 
0.019 and 0.50 ±0.09, respectively. We address the degener¬ 
acy between eccentricity and argument of periapsis following 
Price et al. (2015). Our eccentricities do not vary signifi¬ 
cantly for |w| < 150° away from center-of-transit, consistent 
with Price et al. (2015) and J. Barnes (2007). To find the 
lower limit for eccentricity, we set the center-of-transit at 
periapsis for both planets fit for the eccentricities using our 
gravity-darkened model.we analyze the plausibility of our 
eccentricity values in §5.2. 

4.5. Spin-Orbit Alignment 

Our gravity-darkened model results in a degeneracy in the 
sky-projected alignments between the values A and 180 — A 
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Figure 4. Best-fits and residuals of the KOI-89.01 and KOI-89.02 lightcurves. Red represents the gravity-darkened model, and blue represents 
grazing transits in the spherical model. The KOI-89.02 data are significantly noisier because of KOI-89.02’s much longer orbital period, resulting 
in fewer total transits. We fit the two lightcurves simultaneously, resulting in a single best-fit line. The two lightcurves were placed side-by-side 
for visual comparison. 


(Ahlers et al. 2014). We assume a prograde orbit for KOT 
89.01, constraining A to a single value. This allows us to 
produce single, nondegenerate values for the obliquity (-ip) 
and each planet’s inclination (i) in our best-fit model, which 
we allow to float in the full range of 0° — 360°. 

We find that KOT89 is highly misaligned with a stellar obliq¬ 
uity Ip of 69°±3°, inclinations ii and *2 of 89.340°±0.05° and 
90.64° ±0.06° respectively, and sky-projected alignments Ai 
and A2 of —32° ± 11° and —32° ±40° respectively. The high 
uncertainty of A2 is due to the apparent lack of asymmetry in 
KOl-89.02’s lightcurve because of its photometrically impre¬ 
cise data. We show our constraint of A2 in §5, which removes 
prograde/retrograde degeneracy via dynamic stability tests. 

With these constraints we calculate the true spin-orbit mis¬ 
alignments (fii and ip 2 from the equation (Winn et al. 2007), 

cos((/5i) = sin(^) cos(ii) ± cos{ip) sin(zi) cos(Ai) (1) 

modified for our parameter definitions. We calculate spin- 
orbit alignment angles of 72° ± 3° and 73°ji5^ for the two 
planets respectively. Table 2 lists all of KOT89’s parameter 
constraints. 



KOI-89.02 Transit Time (Hours) 

Figure 5. Lightcurve of both planets transiting simultaneously. Dur¬ 
ing the double transit, the depth is roughly double of a typical transit. 
Our best-fit model (in red) correctly reproduces the time of this event 
and the general shape of the lightcurve it produces. 


4.6. Double Transit 

KOT89.01 and KOT89.02 simultaneously transit halfway 
through 2011, causing a significantly larger transit depth. 
We indicate this event with the arrow in Figure 1 and show 
the double transit and its synthetic lightcurve in Figure 5. 
Our best-fit parameters produce a synthetic lightcurve that 
adequately models this event. We do not find evidence of a 
mutual event in the Kepler dataset. 


Pr, 



( 2 ) 


We find that the star is rotating at 55% — 76% of its breakup 
speed by calculating its rotation period, listed in Table 2: 

p ^ 27ri7„ cos(^) 

V sin(i) 


4.7. Oblateness 

We test the feasibility of KOI-89’s oblateness value of 
0.19^003 by analyzing its breakup rotation period. 


This explains KOI-89’s highly oblate shape and its gravity- 
darkened gradient, which are discussed in §5.3. This star’s 
oblateness is comparable to the fast-rotator Achernar with 
oblateness ~ 0.36 (Carciofi et al. 2008) or other well-known 
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oblate stars such as Altaic (~ 0.2) (Monnier et al. 2007). 
Hence, our model produces physically plausible stellar pa¬ 
rameters. 

5. DYNAMIC STABILITY 

5.1. Coalignment Integrations 

Equation 1 gives a planet’s spin-orbit alignment (pi depen¬ 
dence on the sky-projected alignment A^. We fit for KOI-89’s 
\i in our gravity-darkening model, but are unable to resolve 
A 2 due to its low photometric resolution and the host star’s 
high stellar obliquity ■0. With gravity-darkening-driven 
asymmetry absent in KOI-89.02’s lightcurve, we could not 
fully constrain its transit geometry. 

To estimate KOI-89.02’s sky-projected alignment, we tested 
the system for dynamic stability for various transit geome¬ 
tries. Using our gravity darkening model, we constrained 
KOI-89’s orbital elements with various assumed A 2 values. 
We then used the orbit integrator Mercury from Chambers 
(1999) to test each orbit geometry for dynamic stability. 

Using Mercury, we perform mixed-variable symplectic 
(MVS) integrations of KOI-89 over 10® years using 0.5 day 
timesteps. Using a spherical star allows for physically sound 
integrations that obey the conservation of angular momen¬ 
tum with minimal sacrifice; the stellar J2 ~ 10“^ value (cal¬ 
culated following Murray & Dermont (2008)), coupled with 
the planet’s long orbit periods, cause nodal precession on 
a timescale that would not significantly affect the system’s 
stability. We assume ice-giant densities of p = 1.64 g/cm® 
for both planets. 
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Figure 6 . Survival times for various initial configurations of the KOI- 
89 system. Darker color indicates longer longer survival time, with 
the longest survival time 5.1 x 10^ years. Beyond a. = 20°, ejec¬ 
tion/collision events occur very quickly for all initial configurations, 
suggesting that the KOI-89 system is more stable near coalignment. 



We define an angle a of coalignment between the two orbits, 
defined relative to their angular momentum vectors: 


a = cos 


(i^\ 

vi^iii^2i; 


(4) 


Figure 6 shows the survival time of KOI-89 as a function of 
the coalignment angle a and conjuction longitude. By vary¬ 
ing KOI-89.02’s longitude of periapsis, we vary the conjuc¬ 
tion longitude between the two planets. We define system 
instability as a planet ejection or collision event. None of 
our 360 simulations produced a stable orbit for 10® years, 
indicating we have not found a physically viable system yet. 
In general, survival times are longer for lower a, but the 
longest-lived architectures are non-planar. 

The coalignment angle a is approximately the difference 
between the two planets’ sky-projected alignment angles. 
Using the results of our orbital integrations, we estimate 
the difference between the sky-projected alignment angles 
a ~ IA 2 — All to be 20° ± 20°. This is a conservative esti¬ 
mate based on our results in Figure 6 ; follow-up observations 
would provide a much better calculation of this parameter. 

Our 1332 orbit integrations resulted in a maximum survival 
time of 5.1 x 10^ years. The lack of stable configurations 
suggests that this system is not yet fully understood. If 
the system is in resonance and is non-planar, it may evolve 
chaotically (R. Barnes et al. 2015), and hence long-lived con¬ 
figurations may only exist in small “islands” of parameter 
space. Alternatively, KOI-89’s stability could be brought 
about by unknown additional bodies in the system. We 
show in §5.2 that KOI-89 could be stable if KOI-89.02’s ec¬ 
centricity is lower than our best-fit value of 0.50 ± 0.09. A 
better characterization of this system’s stability could be 
understood via TTV analysis or Rossiter-McLaughlin mea¬ 
surements, but such work is outside the scope of this project. 


5.2. Eccentric Integrations 

In addition to our coalignment/mean longitude stability 
tests, we also test the stability of KOI-89.02’s eccentricity 
of 0.50 ± 0.09 in a coplanar configuration. Van Eylen & 
Albrecht (2015) demonstrated that, in general, multiplanet 
systems have low eccentricities, making KOI-89 a potential 
exception to the rule. See §4.4 for an explanation of our 
treatment of longitude of periapsis. 

We perform a series of integrations in Mercury (Chambers 
1999) using assumed 62 values ranging from 0.0 to 0.95 and 
a 0.05 step size. All 62 < 0.35 are stable, roughly consis¬ 
tent with Petrovich (2015). We show the results of these 
integrations in Figure 7. 

Not surprisingly, lower 62 values yield longer lifetimes and 
overall higher stability. This result suggests three possibili¬ 
ties. The first (and least likely) possibility is that this system 
is in fact not coplanar. If the fitted 62 value of 0.50 ±0.09 is 
correct, then perhaps higher stabilities are found in slightly 
non-coplanar orbits. While higher stability in such a config¬ 
uration is counterintuitive, it does at least reduce the odds 
of a close encounter between the two planets, limiting the 
chances of a violent collision/ejection event. 

The second possibility is that our eccentricity measurement 
contains systematics. A grazing transit would reduce the 
transit duration time similarly to an eccentric orbit transit¬ 
ing near periapsis, and could produce a V-shaped lightcurve 
like we see in Figure 4. KOI-89.02’s low signal-to-noise ratio, 
coupled with the degeneracy between impact parameter and 
planet radius that arises in all grazing transits, prevents us 
from resolving whether KOI-89.02 is in fact fully eclipsing its 
host star. KOI-89.02’s TTVs could also drive up our eccen¬ 
tricity measurement if they are not fully accounted for (Van 
Eylen & Albrecht 2015). High-precision follow-up photom¬ 
etry could better determine KOI-89.02’s orbit parameters. 
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KOI-89.02 eccentricity 

Figure 7. KOI-89 survival times using various assumed 62 values in 
a coplanar configuration. Our N-body integrations are stable through 
10® years for all 62 < 0.35, which is less than two a of our best-fit 
value. The dashed line marks our best-fit value of 0.50 it 0.09. 

including its eccentricity. 

The third possibility is that unknown bodies in the sys¬ 
tem provide stability to these orbits. Antoniadou & Voy- 
atzis (2015) demonstrated that highly eccentric orbits in or 
near mean-motion resonance can exhibit long-term stabil¬ 
ity. Additional bodies could help stabilize KOI-89.01 and 
KOI-89.02, explaining why our best-fit parameters do not 
display dynamic stability through 10® years in our orbit in¬ 
tegrations. 

6 . DISCUSSION 

The KOI-89 system is highly misaligned with spin-orbit 
alignment angles of 72° ± 3° and 73°l'l5^ for the two plan¬ 
ets respectively. Our preliminary dynamical analysis of the 
system §5 failed to find a stable solution, so at this time we 
cannot rigidly constrain the mutual inclination. However, 
we recognize that survival times are longer in general for 
lower a. 

Of ten misalignment mechanisms suggested in the literature, 
five are consistent with our results. We rule out the other 
five mechanisms based on three criteria: 

1. Consistency with KOI-89’s fundamental parameters 
such as orbit period, stellar radius, etc. 

2. The capability to cause extreme misalignment 

3. Conformance with mutually aligned planets 

We compare our results to each mechanism in Table 3. 

6 . 1 . Star-Disk-Binary Interactions 

Batygin (2012) first showed that a stellar companion could 
warp a star’s protoplanetary disk into misalignment. Plan¬ 
ets could then form in the plane of the disk, resulting in 
primordial spin-orbit misalignment (Lai 2014; Xiang-Gruess 
& Papaloizou 2014; Batygin 2012). This mechanism requires 
an unknown binary star in the KOI-89 system, but funda¬ 
mentally agrees with our results in that it could produce 
highly misaligned, coplanar orbits. 


A planet in the potential of a warped protoplanetary disk 
can be driven to very high misalignment values (Terquem 
2013). Teyssandier et al. (2013) found that Jupiter-mass 
planets misaligned from a warped disk experience dynamic 
friction that realigns the planet in timescales shorter than 
the lifetime of the disk. However, Neptune-mass planets can 
remain misaligned and have their eccentricities driven up by 
orbital perturbations from the disk’s gravitational poten¬ 
tial. This mechanism has only been applied to single-planet 
systems, so criterion 3 is inconclusive. However, this mech¬ 
anism agrees with the first two criteria and cannot be ruled 
out based on our results. 

6 . 3 . Planet-Planet Scattering 

Our results cannot entirely rule out planet-planet scattering, 
which is orbit migration due to close encounters between 
high-mass objects (e.g. Chatterjee et al. 2008; Ford et al. 
2005; Raymond et al. 2008; Nagasawa et al. 2008). With 
KOI-89’s net orbital angular momentum highly misaligned 
from the star’s spin angular momentum, conservation of an¬ 
gular momentum would require additional planet (s) to scat¬ 
ter the known two planets. In this scenario it is highly un¬ 
likely that the two planets would end up near mutual align¬ 
ment. However, if a sufficiently large unknown body exists in 
this system, then our orbit integrations are unsound and our 
coalignment constraint for this system is invalid. We there¬ 
fore deem this mechanism consistent as a possible cause of 
KOI-89’s misalignment. Further studies of KOI-89’s TTVs 
could confirm the existence of additional planets. 

6 . 4 . Kozai Resonance 

Kozai resonance in the KOI-89 system requires an unknown 
body that is significantly misaligned with its known orbital 
plane (Libert & Tsiganis 2009; Thies et al. 2011; Payne et al. 
2010). Such an event would likely not produce coplanar 
orbits for KOI-89.01 and KOI-89.02. However, Kaib et al. 
(2011) suggests that coplanar, inclined orbits might arise as 
a result of this mechanism. We therefore deem this method 
consistent. 

6 . 5 . Internal Gravity Waves 

Rogers et al. (2012) showed that angular momentum trans¬ 
port between the convective interior and radiative exterior 
of hot, early-type stars can change the observed stellar spin 
axis, resulting in spin-orbit misalignment. This misalign¬ 
ment mechanism happens independently of orbiting bodies 
and does not affect coplanarity. The 2-D simulations per¬ 
formed in Rogers et al. (2012) found that this mechanism 
can occur on a timescale as short as tens of years, and can 
explain retrograde orbits. Whether this mechanism can pro¬ 
duce spin-orbit misalignments near 90° is still under inves¬ 
tigation. 

6 . 6 . Planet-Embryo Collisions 

Planet-embryo collisions can occur in any standard forma¬ 
tion model, and they can drive migration in various ways 
(Levison et al. 1998; Charnoz et al. 2001). However, this 
mechanism can produce large spin-orbit misalignment angles 
only for small rocky bodies and does not apply to the highly 
misaligned giant planets KOI-89.01 and KOI-89.02. Addi¬ 
tionally, this mechanism likely could not produce coplanar 
misaligned orbits because the collisions driving this mecha¬ 
nism are unique to each planet. 


6 . 2 . Inclination Driven By A Warped Disk 


6 . 7 . Chaotic Evolution of Stellar Spin 




Mechanism 

Relevant Works 

Criteron 1 

Criteron 2 

Criteron 3 

Star-disk-binary interactions 

Lai (2014), Xiang-Gruess & Papaloizou (2014), 
Batygin (2012) 

Yes 

Yes 

Inconclusive 

Inclination driven by a warped disk 

Terquem (2013) 

Yes 

Yes 

Inconclusive 

Planet-planet scattering 

Chatterjee et al. (2008), Ford et al. (2005), Ray¬ 
mond et al. (2008), Nagasawa et al. (2008) 

Yes 

Yes 

Inconclusive 

Kozai resonance 

Libert &: Tsiganis (2009), Thies et al. (2011), 
Payne et al. (2010) Kaib et al. (2011) 

Yes 

Yes 

Inconclusive 

Internal Gravity Waves 

Rogers et al. (2012), Fuller et al. (2015) 

Yes 

Yes 

Inconclusive 

Planet-embryo collisions 

Levison et al. (1998), Charnoz et al. (2001) 

Yes 

No 

No 

Chaotic evolution of stellar spin 

Storch et al. (2014), Valsecchi &: Rasio (2014) 

No 

Yes 

Yes 

Magnetic torquing 

Lai et al. (2011), Spalding & Batygin (2014) 

No 

Yes 

Yes 

Coplanar high-eccentricity migration 

Petrovich (2014) 

Yes 

No 

Yes 

Inclination Resonance 

R. Barnes et al. (2015) 

No 

Yes 

Inconclusive 


Table 3 

Possible spin-orbit misalignment mechanisms for the KOI-89 system. We list ten possible causes of spin-orbit misalignment that have been put 
forward in the literature and rule out five of them based on our best-fit results and our estimation of the two planets’ coalignment. The three 
criteria are: (1) consistency with KOI-89’s fundamental parameters, (2) the capability to cause extreme misalignment, and (3) the production of 

mutually aligned planets. 


Storch et al. (2014) and Valsecchi & Rasio (2014) demon¬ 
strated that strong tidal dissipation can cause chaotic evo¬ 
lution of stellar spin. This mechanism requires hot Jupiters 
with periods < 3 days. Such a body in the KOI-89 system 
would have to be drastically misaligned from the plane of the 
other two orbits; therefore, this mechanism cannot be the 
standalone cause of misalignment because some other mech¬ 
anism would have to misalign the hot Jupiter. If there was 
a non-transiting hot Jupiter that was initially misaligned, 
it could torque the star into misalignment with the other 
planets. If said hot Jupiter fell into its host star because 
of tidal decay, it could change both KOI-89’s rotation axis 
and rotation rate (Jackson et al. 2009). However, early-type 
stars such as KOI-89 have weak tidal interactions in general 
(Ogilvie & Lin 2007). 

6 .8. Magnetic Torquing 

Magnetic torquing between a stellar magnetic field and a 
protoplanetary disk can cause misalignment by torquing the 
disk away from the star’s equatorial plane (Lai et al. 2011; 
Spalding & Batygin 2014). KOI-89 is an early-type star 
with a weak magnetic field (Bagnulo et al. 2002), so this 
mechanism could not cause KOI-89’s high misalignment. We 
note that the magnetic fields of fast-rotators are still under 
investigation (Ibahez-Mejia & Braithwaite 2015); a better 
understanding of these magnetic fields may reveal this to be 
a possible misalignment mechanism for KOI-89. 

6.9. Coplanar High-Eccentricity Migration 

Coplanar high-eccentricity migration can occur in mutually 
aligned multiplanet systems with at least one highly ec¬ 
centric orbit. Secular gravitational effects excite the inner 
planet’s eccentricity to very high values, and planetary tidal 
dissipation during periapsis reduces the orbit’s semi-major 
axis. This mechanism occurs primarily in the planets’ or¬ 
bital plane, predominantly maintaining the system’s original 
spin-orbit alignment angles (Petrovich 2014). 

6 .10. Inclination Resonance 

If the planets are in resonance and possess a mutual inclina¬ 
tion, then the orbital inclinations can be driven to very large 
values (R. Barnes et al. 2015). In that case we may expect 
to find at least one planet in a misaligned orbit. This phe¬ 
nomenon can also produce very large eccentricities. How¬ 
ever, this mechanism depends on stellar torquing from tidal 
interactions for both planets to be discovered in a misaligned 


state. Such tidal interaction is weak around early-type stars 
(Ogilvie & Lin 2007). While this is a possible cause of KOI- 
89’s extreme misalignment, it requires an external mech¬ 
anism to bring about an initial mutual inclination. More 
work is needed to understand if this scenario is possible and 
could apply to KOI-89. 

7. CONCLUSION 

We constrain the individual spin-orbit alignments of multi¬ 
planet system KOI-89. With our gravity-darkened model, 
we found significant spin-orbit misalignment with angles of 
72° ± 3° and 73°^5“ for KOI-89.0I and KOI-89.02, respec¬ 
tively. We also constrain other fundamental parameters of 
the KOI-89 system and estimate the mutual alignment be¬ 
tween KOI-89.01 and KOI-89.02. We show these results in 
Table 2. 

While our measurements alone do not uniquely assign a mis¬ 
alignment mechanism to KOI-89, the large spin-orbit align¬ 
ment angles ipi and low coalignment angle a of this system 
limit the possible mechanisms for planet migration. These 
values, the measured TTVs, the near 5:2 resonance, and the 
fast rotation of the star itself all imply a dynamic formation 
history. 

KOI-89 is of of particular interest because it can experimen¬ 
tally constrain the numerous outstanding hypotheses that 
have been proposed to generate misalignment. We limit 
possible causes to star-disk-binary interactions, disk warp¬ 
ing via planet-disk interactions, planet-planet scattering, or 
internal gravity waves in the convective interior of the star. 
Follow-up observations searching for additional bodies could 
provide evidence for any of these hypotheses, including in¬ 
ternal gravity waves if no additional bodies are found. 

Much could still be learned about the KOI-89 system. As- 
teroseismic determination of the star’s oscillation modes 
could confirm various stellar properties such as stellar ra¬ 
dius, mass, and obliquity. Doppler tomographic observa¬ 
tions could constrain the individual ascending nodes of the 
two planets. Analysis of the TTV could confirm/deny the 
existence of undiscovered planets in the system. High- 
precision photometry could better constrain the two planets’ 
eccentricities and impact parameters and help resolve the 
degeneracy between limb-darkening and gravity darkening. 
Any of these follow-up observations would shed new light on 
the formation of solar systems dissimilar to our own. 

The constraints provided in this work add to the sample of 
known misaligned systems - particularly misaligned multi- 
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planet systems, of which only a small nnmber are currently 
known. The unique nature of the KOI-89 system provides 
new insight for studying system formation and evolution. It 
also adds to the surprising diversity of exosystems discov¬ 
ered to date. Future studies can apply the knowledge gained 
from this work to a wide variety of misaligned and dynamic 
systems. 

We would like to thank the Kepler Science Team for mak¬ 
ing this work possible - particularly Dr. David Latham for 
providing the nsin(z) measurement of this system. JPA and 
JWB are funded by NASA Proposal 7^13-ADAP13-213. RB 
acknowledges support from NSF grant AST-1108882. 
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